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ABSTRACT 



We report the detection of broad Ho- emission in three X-ray selected obscured AGNs at z~l-2. By exploiting the Her width and 
the intrinsic X-ray luminosity, we estimate their black hole masses, which are in the range 0.1-3xl0 9 M Q . By means of multi-band 
photometric data, we measure the stellar mass of their host galaxy and, therefore, infer their M BH /M star ratio. These are the first 
obscured AGNs at high-z, selected based on their black hole accretion (i.e. on the basis of their X-ray luminosity), that can be located 
on the M BH -M star relation at high-z. All of these obscured high-z AGNs are fully consistent with the local M BH -M slar relation. This 
result conflicts with those for other samples of AGNs in the same redshift range, whose MBH/M slar ratio departs significantly from the 
value observed in local galaxies. We suggest that the obscured AGNs in our sample are in an advanced evolutionary stage, have already 
settled onto the local MBH-M star relation, and whose nuclear activity has been temporarily revived by recent galaxy interactions. 
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1. Introduction 

A major breakthrough in our understanding of galaxy evolution 
has been the discovery of a tight correlation, in the local uni- 
verse, between the mass of supermassive black holes (Mbh) and 
the mass of their host spheroids (e.g. Magorrian et al. 1998, 
Ferrarese & Merritt 2000). The existence of this relation im- 
plies a strong physical connection between galaxy formation and 
growth of black holes at their centers. Various models and sim- 
ulations have been proposed to explain this correlation (e.g., 
Menci et al. 2006, Marulli et al. 2008, Hopkins et al. 2006, 
Volonteri & Natarajan 2009). These models predict different 
evolutionary patterns on the Mbh _ M star diagram through the 
cosmic epochs. Therefore, determining the location of high-z 
galaxies on the Mbh - M stal - diagram provides a crucial test for 
galaxy-black hole coevolutionary scenarios. 

Various observational studies have attempted to measure the 
evolution of the Mbh - M stal - relation at high redshift. Most of 
these studies infer the Mbh in high-z AGNs by using "virial 
estimators", involving the AGN luminosity and width of the 
"broad emission lines" (e.g. Vestergaard et al. 2009, Shen et al. 
2010). Broad lines, especially in the UV (observed in the opti- 
cal at high-z), are generally detected only in unobscured, type 
1 AGNs (AGN Is), whose broad line region can be observed di- 
rectly. Thus the investigation of the Mbh - M stal - relation has 
focused mostly on AGN Is, especially at high-z. Jahnke et al. 
(2009) use the virial relations to measure the BH masses of ten 
type 1 AGN at z~1.4 and derive the stellar masses of their host 
galaxy by means of multi band fitting. They find that the BH- 
to-total stellar mass ratio does not evolve relative to the local 
relation; however, since several hosts show evidence of a sub- 



stantial disk component they suggest that the ratio of BH mass to 
bulge mass (given by the difference between total and disk stellar 
mass) probably evolves with redshift. Merloni et al. (2010) ex- 
pand the Jahnke et al. sample to 89 X-ray selected type 1 AGN at 
l<z<2.4 in COSMOS (adopting a different method to measure 
the stellar masses). They infer that the MBH/M star ratio increases 
significantly, by a factor of about two, relative to the local rela- 
tion. A similar result was obtained by Peng et al. (2006) by using 
a sample of lensed quasars in a similar redshift range. At higher 
redshifts (z~4-6), dynamical mass measurements, obtained by 
exploiting CO maps, suggest that the MBH/M star ratio in type 1 
AGNs increases even further, up to an order of magnitude rel- 
ative to the local relation (Walter et al. 2004; Maiolino et al. 
2007a, Lamastra et al. 2010, Wang et al. 2010). 

A possible caveat of these observational studies is that the 
targets are extracted from AGN-selected samples, i.e. on the ba- 
sis of the black-hole accretion rate. This is expected to introduce 
a bias in favor of massive black holes (which can reach higher 
absolute accretion rates within their Eddington limit) Q More 
specifically, even if the Mbh _ M star relation does not evolve with 
redshift, it does have a scatter and the AGN-selection bias fa- 
vors the selection of objects with MBH/M stal higher than the true 
distribution, hence mimicking an evolution. Lauer et al. (2007) 
estimate that this bias may increase the inferred Mbh _ M star 
even by a factor of a few, depending on the AGN luminosity and 
the intrinsic scatter of the Mbh - M stal relation at high redshift. 
However, Merloni et al. (2010) demonstrate that, at least for their 
X-ray selected, unobscured AGN 1 sample, this bias should not 
significantly affect their own results. 



* Based on data obtained at the VLT through the ESO programs 73 A- 
0598, 076.A-0681, 077.B-0368. 



1 Quasars with dynamical masses inferred from CO maps are also 
selected based on their millimeter continuum, hence for this subsample 
the bias may be more complex. 
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Fig. 1. Near-IR spectra (zoomed around Ha) of the three X-ray selected obscured AGNs with evidence for a broad component 
of Ha. The red line shows the fit resulting from a combination of a broad Ho- Gaussian component (blue line), narrow Gaussian 
components of Ha, [Nil] and, in the case of Abell2690#29, [SII] (green lines) and a linear continuum. The bottom panels show the 
sky spectrum in arbitrary units. 



Table 1. Properties of the X-ray selected obscured AGN with broad Ha. 



ID 


z 


instr. 


logl^-lOkev" 


logN H 


logLboi" 


FWMVW' 


logMflH 


L/LEdd 


logM star 








[erg s~'] 


[crrr 2 ] 


[erg s~'] 


[kms- 1 ] 


[M ] 




[M Q ] 


Ab2690#29 


2.087 


ISAAC 


44.95 


22.80 


46.81 


5871±245 


9.44±0.31 


0.12 


12.50±0.20 


XMMES1_460 


1.748 


SINFONI 


44.86 


22.50 


46.59 


2316+113 


8.54±0.31 


0.63 


11.53±0.20 


BPM1627#181 


1.335 


SINFONI 


44.20 


22.81 


45.76 


2491+511 


8.15±0.36 


0.32 


11.28+0.20 



Notes. <o) Absorption corrected 2-10 keV luminosity. (W Bolometric luminosity inferred by using the X-ray bolometric corrections given in 
Marconi et al. (2004). (c) Full Width at Half Maximum of the broad component of Ha in km/s. 



The only work investigating the evolution of the Mbh - M stal - 
relation through a sample of objects not pre-selected among 
AGNs is that of Alexander et al. (2008). Their parent sam- 
ple is composed of sub-millimeter galaxies (SMGs), which 
are strongly starbursting systems at z~2. Their sample consists 
mainly of obscured AGNs hosted such systems. For a subsam- 
ple of these obscured AGNs, which show some broad Ha or H/3, 
they are able infer the BH mass. In contrast to the type 1 AGN 
studies mentioned above, the inferred MBH/M star ratio of SMGs 
at z~2 is lower than observed locally. However, this sample may 
be subject to an opposite bias than AGN-selected samples. Given 
the extraordinarily high star-formation rate of SMGs and the 
correlation between SFR and stellar mass, the sample used by 
Alexander et al. (2008) may be biased toward objects of high 
stellar mass (although SMGs are outliers on the SFR-mass rela- 
tion, toward higher SFRs). The sub-mm selection may also favor 
dust-rich systems (Santini et al. 2010), hence more-evolved host 
galaxies and more massive galactic hosts. Altogether, these ef- 
fects may favor of objects with MBH/M stal - ratios lower than the 
real distribution. 

Here we present near-IR spectra of three obscured quasars at 
z~l-2, extracted from a sample of hard X-ray selected AGNs, 
that display a broad component of Ha, allowing us to infer 
the BH mass and MBH/M star ratio. These are the first obscured 
AGNs at high-z selected on the basis of their black hole ac- 
cretion (i.e. because X-ray selected) that can be located on the 
Mbh - M stal diagram at high-z. We illustrate that the trend of 
these obscured AGNs differs significantly from that of the other 
samples of AGNs investigated in the same redshift range. 

2. Sample selection, observations, and data 
analysis 

These near-IR spectroscopic observations are part of a program 
designed to determine the redshift of and characterize obscured 
AGNs at high redshift selected from wide-area hard X-ray sur- 



veys. These AGNs are too faint at optical wavelengths to be ob- 
served spectroscopically. The parent samples consist of XMM 
hard X-ray sources from the HELLAS2XMM extended survey 
(Cocchia et al. 2007) and the XMM survey of the ELAIS-S1 
field (Feruglio et al. 2008). From these samples, we selected 14 
sources characterized by a X-ray to optical flux ratio X/O>10, 
which is typical of most obscured (Compton-thin) AGNs at high 
redshift and, in particular, obscured QSOs (Fiore et al. 2003). 
Most of these optically faint AGNs have also very red optical- 
to-near-IR colors and, more specifically, R - K > 5 (Mignoli 
et al. 2004), i.e. are extremely red objects (EROs). This class 
of objects is generally found to consist mostly of either quies- 
cent, evolved galaxies or dust-reddened star-forming galaxies at 
z~l-2, with a fraction of obscured AGNs (e.g. Campisi et al. 
2009). In the specific case of X-ray selected EROs, Mignoli et 
al. (2004) find that most of them are resolved in the K-band with 
typical sizes of 0".5 and elliptical-like profiles, implying that the 
K-band light is generally dominated by an early type host. The 
latter result is also confirmed in most objects by a SED analysis 
(Pozzi et al. 2007, 2010). 

Our observations were performed in three observing runs. 
In the first run, we used ISAAC at the VLT to observe three 
HELLAS2XMM sources. We used ISAAC in its low spectral 
resolution mode (R~500), covering the bands J, H, and K bands. 
These observations were presented in Maiolino et al. (2006) and 
we refer to that paper for further details. In the second run, we 
used SINFONI, the near-IR integral field spectrometer at the 
VLT, to observe three HELLAS2XMM sources. We used both 
the H+K R-1500 grating and the J-band R-2000 grating. The 
on-source integration times ranged from 40 min to 4 h per target 
in each band. In the third run, we observed nine sources from the 
ELIAS-S1 sample with SINFONI, with the same instrumental 
setting as for the previous SINFONI run and similar integration 
times. The seeing during the observations ranged from 0".7 to 
1".5. Data reduction was performed using the ESO pipeline. 
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Fig. 2. M B h - M star relation for AGNs at z~l-2. The red dia- 
monds with errorbars show the location of the three X-ray se- 
lected obscured AGNs presented in this paper. Green triangles 
are type 1, X-ray selected unobscured AGNs in COSMOS from 
Merloni et al. (2010). Blue squares are obscured AGNs in SMGs 
from Alexander et al. (2008) with broad Ha and stellar masses 
from Hainline et al. (2010). The black line is the local relation 
obtained by Sani et al. (2010). 



All targets were detected in the continuum. However, only 
for seven of them could we identify emission lines that al- 
low us to unambiguously determine their redshift. At the in- 
ferred redshift (l<z<2.1) the inferred X-ray luminosity and X- 
ray spectral shape imply that these are all type 2, obscured QSOs 
(L2_iokeV > 10 44 erg s -1 , Nh > 10 22 crrT 2 ). Here we focus on 
three sources for which we detect a broad component of Ho- that 
allowed us to infer their BH masses. Their main observational 
and physical properties are reported in Table 1 . 

3. The M B h - M star relation 

In Fig. Q~| we show the near-IR spectra of the three obscured 
quasars at z~l-2 with broad Ha (zoomed around the Ho- line). 
No H/3 broad line is detected for any of these sources, which 
is indicative of a type 1.9 optical classification and consistent 
with the absorption inferred from both their X-ray spectrum and 
their optical faintness. We fitted the line profile with a broad 
Gaussian and three narrow Gaussians (FWHM<800 km/s) ac- 
counting for the contribution of the narrow components of both 
Ha and [NII]6548,6584. In the case of Abell2690#29, we also 
included the [SII]6720 doublet because the Ha is so broad that it 
extends to this rest-frame wavelength. We assumed that the nar- 
row components have the same width and velocity shift. We also 
assumed that F([NII]6584)/F([NII]6548)=3, as given by the rela- 
tive transition probabilities of the corresponding levels. The con- 
tinuum in this narrow redshift range is fitted simply with a linear 
component. Regions affected by very low S/N or bad residuals 
due to strong sky lines or bad atmospheric transmission were ex- 
cluded from the fit. The resulting spectral fits are shown in Fig.Q] 



and the width of the broad components are listed in Table 1 . The 
uncertainty in the width of the broad line is also estimated by 
accounting for the correlation with the other components. 

We cannot use the virial formula for the BH mass based on 
the Ha FWHM and luminosity proposed by Greene & Ho (2005) 
since the Ha is likely absorbed, hence its intrinsic luminosity is 
difficult to measure. However, we can use both the virial for- 
mula for the BH mass involving the optical continuum luminos- 
ity AL A (51Q0A) and FWHM^s (Marconi et al. 2010) 



M BH = 6.16 x 10 6 



^(5100) 
10 44 ergs-' 



0.5 



FHWMh/; 
10 3 km s- 1 



MfT 



(1) 



and the relation between L2_iokeV and ^IL^ (5 1 00) (Eq.5 in 
Maiolino et al. 2007b) to estimate the optical continuum, and 
replace FWHM^ with FWHMh» by using the relation given in 
Greene & Ho (2005) (Eq. 3 therein). By replacing these relations 
in Eq.[T| we obtain 



M BH = 1.56X10 7 



l-2-10keV 

10 44 erg s" 



0.694 



FHWMho 
10 3 kms- 1 



M Q .(2) 



Since we know the absorption-corrected 2-10 keV luminosity of 
our three sources (Table 1), the width of the broad Ha inferred 
from our near-IR spectra allows us to infer the BH by exploiting 
Eq.|2] The inferred BH masses are reported in Table 1, and range 
between 1.3 x 10 8 M and 2.8 x 10 9 M . The uncertainty in 
Mbh is dominated by the dispersion in the Mbh - c, Lt, u i local 
relation (0.3 dex), which is used to calibrate the virial relations. 
The uncertainty in the L2-iokeV - L5100 relation coefficients also 
contributes with 0.07 dex (in quadrature), while the uncertainty 
in the width of the broad lines contributes significantly only in 
BPM1627#181 with 0.19 dex (in quadrature). 

For consistency, we checked that the masses inferred by us- 
ing the Ha luminosity and the relation proposed by Greene & 
Ho (2005) imply lower BH masses, confirming that Ha is ab- 
sorbed. Similarly, the direct use of the relation in EqQ]involving 
L / |(5100), by taking the observed optical (B rest-frame) contin- 
uum of our sources, deduces BH masses lower than those pre- 
sented in Table 1, confirming that the AGN optical light is ab- 
sorbed and that the observed rest-frame blue continuum is dom- 
inated by the host galaxy (as inferred from the SED fitting). 

The BH masses measured by ourselves imply that these three 
obscured AGNs are accreting at a rate that is about 0.1 - 0.6 of 
their Eddington limit (Table 1). 

The stellar masses were inferred by exploiting the multi- 
wavelength photometric points available for these sources and 
the SED fitting code (Fritz et al. 2006) used in Pozzi et al. (2010) 
and Vignali et al. (2009) for a subsample of obscured AGN in 
the HELLAS2XMM sample. This code combines synthetic stel- 
lar libraries, AGN dusty torii models (Fritz et al. 2006) including 
the intrinsic AGN component (if visible), and starburst IR tem- 
plates to reproduce the observed SED. We refer to Pozzi et al. 
(2010) for a detailed description of the method. Here we only 
emphasize that in the case of obscured AGNs, the rest-frame op- 
tical to near-IR radiation is generally dominated by the stellar 
light from the host galaxy and therefore constraining the stel- 
lar mass is easier than for the host galaxies of type 1 AGNs. The 
stellar masses inferred for the host galaxies of the three obscured 
quasars presented here are reported in Table 1 . We note that to 
compare with the local relation, which is between the BH mass 
and the stellar mass of the spheroid, we would need to extract 
the spheroidal component of the stellar mass in our high-z tar- 
gets. However, Mignoli et al. (2006) demonstrate that the class 
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of targets in our sample is generally characterized by elliptical- 
like profiles, hence the bulk of the stellar light is associated with 
a massive spheroid. 

The resulting location of the three X-ray selected, obscured 
AGNs at z~l-2 is shown in Fig. [2] (red diamonds) along with 
the location of type 1, unobscured AGNs in the same redshift 
range from Merloni et al. (2010) (green triangles) and the ob- 
scured AGNs in SMGs at z~2 (blue squares). For the SMG sam- 
ple we only plot those AGNs with broad Balmer lines detec- 
tion (Alexander et al. 2008) and with available stellar masses 
from Hainline et al. (2010). The solid line shows the local re- 
lation deduced by Sani et al. (2010). All of the X-ray selected 
obscured AGNs are fully consistent with the local Mbh _ M star 
relation and do not show any evidence of evolution, in contrast 
to both type 1, unobscured X-ray selected sample and the ob- 
scured AGNs selected among SMGs. For our obscured AGNs 
we find that (log (M B H/M stal )) = -3.06 + 19, which is fully con- 
sistent with the local ratio log (MBH/M star )i ocal = -3.01 obtained 
by Sani et al. (2010) in a similar mass range. 

4. Discussion 

Our sample is small and does not allow us to draw firm conclu- 
sions from a statistical point of view. However, it is very intrigu- 
ing that, in contrast to other high-z AGN samples, the three X- 
ray selected obscured AGNs do not show any systemic deviation 
from the local Mbh - M star relation. If this trend were confirmed 
for a larger sample of X-ray selected obscured AGNs, it would 
have important implications for understand either biases affect- 
ing current studies or the evolution of the MBH/M stal ratio during 
different evolutionary stages of galaxies. 

SMGs have a large scatter, but on average they tend to be be- 
low the local Mbh - M star relation (Alexander et al. 2008). The 
difference with our results may be explained in terms of the bias 
discussed in the introduction, i.e. SMG AGN hosts are biased to- 
ward higher stellar masses because of the mass-SFR relationship 
and/or because of their higher dust content (hence more evolved 
hosts). Alternatively, the two samples may trace different evolu- 
tionary stages. According to many theoretical models, black hole 
accretion occurs predominantly in the phases of galaxy merg- 
ing and interactions (which destabilize the gas, allowing it to 
flow towards the nucleus). Lamastra et al. (2010) suggest that 
SMGs may be objects where previous secular star formation has 
increased the stellar mass, while the lack of previous mergers 
and interactions has prevented the BH mass from growing sig- 
nificantly. This may explain their lower BH-to-stellar mass ra- 
tio. According to this scenario, SMGs are being observed dur- 
ing their first major-merging event. In contrast, X-ray selected, 
obscured AGNs (at least in our ERO subsample) may be the de- 
scendants of SMGs, observed at a later stage when, at the end 
of the whole interaction/merging phase, galaxies have already 
settled onto the local MBH/M stal relation. Our targets have in 
general elliptical-like profiles (Mignoli et al. 2005) and an SED 
typical of quiescent galaxies (Pozzi et al. 2010), therefore sup- 
porting the latter scenario. However, we that a few objects in the 
parent sample have strong far-IR and submm emission indica- 
tive of vigorous star-forming activity (Vignali et al. 2009, Pozzi 
et al. 2010). In the specific case of our three targets, the optical 
to near-IR photometric points do not allow us to provide tight 
constraints on the SFR. 

The differences between the location on the Mbh - M star di- 
agram of high-z X-ray selected obscured AGNs and X-ray se- 
lected unobscured AGNs is more puzzling. Except for our sam- 
ple being obscured, our selection criteria do not differ signifi- 



cantly from the Merloni et al. (2010) sample (e.g. similar X-ray 
luminosity range), hence should not produce differential selec- 
tion effects (in particular, our sample should also be unaffected 
by the Lauer et al. 2007 bias). However, we note that the distri- 
bution of unobscured AGNs has a large scatter. The two high-z 
obscured AGNs with low stellar masses (i.e. overlapping with 
the Merloni et al. sample in terms of stellar mass) may at least 
still be consistent with the tail of the distribution of unobscured 
AGNs in the same redshift range. Larger number statistics are 
required for the high-z obscured sample to clarify whether they 
are, on average, offset or consistent with the unobscured samples 
in terms of black hole-to-galaxy mass ratio. Were future observa- 
tions to confirm the difference between unobscured and obscured 
(X-ray selected) AGNs, this may suggest that the two classes 
are associated with different evolutionary stages. The general 
expectation of models is that unobscured, type 1 AGNs repre- 
sent a later evolutionary stage than to the obscured growth phase. 
However, the finding that unobscured AGN are offset from the 
local Mbh - M star relation and are moving towards it (Merloni 
et al. 2010), while obscured AGN have already settled onto the 
local relation, suggests that the latter are in a later evolution- 
ary stage. A possibility is that the obscured AGNs in our parent 
sample of EROs are in a late phase where the BH and their host 
galaxy have already reached the local Mbh - M stal relation, but 
are temporarily rejuvenated by a late accretion episode, making 
them detectable as quasars. In support of this scenario, we men- 
tion that our data reveal that more than 25% of the targets in our 
parent sample (and all of the sources presented here) are in close 
interaction with one or more galaxies located within a radius of 
~ 15 kpc, which may be responsible for triggering BH accretion. 

A larger sample of high-z obscured, X-ray selected AGNs 
with broad Ha (to infer the BH mass) is certainly required to 
confirm with higher statistical significance that this class of tar- 
gets follow the local Mbh - M star relation and that they differ 
from other high-z AGN samples, as well as to test the various 
possible scenarios discussed above. 
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